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Mirometeoroid Complex: The interaction of the micrometeoroid complex with the lunar 
surface is evidenced by numerous glass-lined microcraters on virtually every lunar surface 
exposed to space. Such craters range in size from < .l/xm to approximately 2 cm in 
diameter. With the use of small-scale laboratory cratering experiments for "calibration," 
the observed crater-sized frequency distributions may be converted into micrometeoroid 
mass distributions. These "lunar” mass distributions are in essential agreement with 
satellite data for masses > 10~“g. However, for masses < 10'“ g there is considerable 
discrepancy. A radiation pressure cutoff does not exist because masses as small as 10'“ g 
can be observed. The absolute flux of micrometeoroids based on lunar rock analyses 
averaged over the past few 10" years is approximately an order of magnitude lower than 
presentday satellite fluxes; however, there is indication that the flux increased in the past 
10* years to become compatible with the satellite data. Furthermore, there is detailed evi- 
dence that the micrometeoroid complex existed throughout geologic time. 

Some physical properties of micrometeoroids may be deduced by comparing lunar crater 
geometries with those obtained in laboratory experiments. The preponderance of circular 
outlines of lunar microcraters necessitates equidimensional, if not spherical, micro- 
meteoroids. Irregular shapes such as whiskers, needles, platelets, rods, etc. — postulated 
in the past — do not contribute' substantially to the micrometeoroid population and are 
rare, if not absent. The depth/diameter ratios of lunar microcraters are compatible with 
micrometeoroid densities of 2 to 4 g/cm’; densities < 1 g/cm" can be excluded. These 
findings have astronomical significance with respect to comets, i.e., the source area for 
micrometeoroids. 

Regolith-Dynamics: Monte Carlo-based computer calculations, as well as analytical ap- 
proaches utilizing probabilistic arguments, were applied to gain insight into the principal 
regolith impact processes and their resulting kinetics. Craters 10 to 1500 m in diameter 
are largely responsible for the overall gowth of the regolith. As a consequence the regolith 
has to be envisioned as a complex sequence of discrete ejecta blankets. Such blankets 
constitute first-order discontinuities in the evolving debris layer. The micrometeoroid com- 
plex then operates intensely on these fresh ejecta blankets and accomplishes some degree 
of mixing and homogenization. True mixing, however, can be accomplished only in an 
uppermost layer of approximately 1-mm thickness, before a new ejecta event covers this 
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layer and effectively removes it from the zone of active reworking. While, e.g., a layer 1 cm 
in depth is turned over only one time in approximately 10' years, the uppermost 1 mm of 
that surface has been turned over already 250 times and the uppermost .1 mm more than 
2000 times during the same period. Therefore the lunar regolith becomes rapidly quiescent 
with depth. Though the micrometeoroid bombardment is extensive, a startigraphic sequence 
may readily be preserved, as evidenced in returned core tube materials. The erosion of 
lunar rocks caused by micrometeoroids is calculated at .3 to .6 mm per 10" years. The 
mean surface residence time of a rock of 1 kg in mass is in the order of 3 x 10" years, be- 
fore it will be catastrophically destroyed by rupturing due to the impact of large micro- 
meteoroids. This catastrophic destruction is far more effective than single particle abrasion 
in obliterating a lunar rock specimen. Due to the vagarities of the random impact process, 
caution is necessary to delineate regolith dynamics from lunar sample analyses that are 
not based on a statistically significant number of observations. 


With increasing resolution of lunar surface 
photographs prior to actual sample return, 
it became more and more obvious that mete- 
oroid impact had played a substantial role in 
the evolution of the lunar surface. It was dis- 
covered that meteoroid impact had operated 
on scales from hundreds of km down to a few 
cm (ref. 1). However, immediately upon cur- 
sory inspection of returned rocks it was 
learned that impact processes also occurred 
on still smaller scales : the ubiquitous presence 
of glass-lined lunar microcraters was ample 
evidence that virtually every lunar surface 
exposed to space was also subjected to the 
bombardment of micrometeoroids. In the 
meantime numerous laboratory investigations 
revealed that many properties of the lunar 
regolith are either directly or indirectly dom- 
inated by impact processes far beyond the 
original expectations. A proper understand- 
ing of many regolith processes therefore de- 
pends critically upon an understanding of 
the regolith impact history. 

A thorough understanding of this history 
is possible only by combining lunar observa- 
tional data, laboratory impact experiments, 
and theoretical calculations. This report at- 
tempts to summarize such analyses. We will 
first discuss observational evidence of lunar 
microcraters and its implications to the 
micrometoroid complex, including some astro- 
nomical, consequences. We then will present 
some analytical and computer-based calcula- 
tions that will aid in the understanding of 
some principal regolith processes as well as 
their kinetics. Due to limited space some 
detailed argumentation cannot be presented. 


and the reader must be referred to the orig- 
inal reports. In addition, a multitude of other 
interesting observations and interpretations 
had to be deleted. We attempted, however, to 
present the most important aspects of the 
impact process as we understand them today. 

Lunar Data of the Micrometeoroid 
Complex 

MICROCRATER-MORPHOLOGY 

Glass surfaces are by far the most suitable 
materials to study micrometeoroid impacts, 
because in comparison with crystalline rocks 
and breccias, they are usually smooth and 
observational conditions are optimized (fig. 
1). Furthermore, glasses are also the best 
investigated materials in small-scale labora- 
tory cratering simulations. Thus — unless 
specified — the detailed morphology data, 
crater-size frequency distributions, and as- 
sociated flux considerations are derived from 
lunar glass surfaces only. 

Microcraters on lunar glass surfaces may 
range in diameter from less than .l/^m up to 
approximately 2 mm; on crystalline rocks 
craters as large as 20 mm in pit diameter 
were observed. Crater morphology differs 
characteristically as a function of absolute 
crater diameter (refs. 2-5). Craters smaller 
than l/xm are cup-shaped, glass-lined depres- 
sions — termed “pit” — with a pronounced rim 
of molten target material (fig. 2(a)). Craters 
between l;um and 10/un in pit diameter (figs. 
2(b) and 2(c) ) are transitional between the 
above morphology and that typical for craters 
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Figure 1 . — Large glass-coating on lunar rock 6UU55 
with abundant microcraters. All structures are 
above 5 jum in diameter and therefore display char- 
acteristic spall zones. Close to the fracture zone 
exposing the underlying anorthositic substrate, the 
crater densities are very high and in approximate 
equilibrium. (Sidelength of picture: 3.2 cm.) 

larger than lO^am. When above lOjum in 
diameter, they not only possess a central 
glass-lined pit but also a concentric spall 
zone (fig. 2(d)). The spall zone may or may 
not be totally spalled off for craters between 
10/im and 50/im, but all craters above SO/un 
in diameter have a completely developed 
spall zone. Morrison et al. (ref. 5) delineated 
the following relationship: Ds — 2.37 X 
where Ds is the spall zone and Dp the 
pit diameter. 

For comparison, identical structures pro- 
duced in the laboratory are illustrated in 
figures 2(e) and 2(f). Laboratory crater 
studies performed by electrostatic particle 
accelerators (refs. 6-13) indicate that a 
glass-lined pit is produced only at projectile 
velocities exceeding 3 km/s. The development 
and extent of a spall zone characteristic for 
the larger lunar craters requires velocities in 
excess of 5 km/s. 

PHYSICAL PROPERTIES OF 
MICROMETEOROIDS 

Laboratory simulations by Mandeville and 
Vedder (ref. 9) ; Kerridge and Vedder (ref. 


14) ; Vedder and Mandeville (ref. 12) ; and 
Mandeville (ref. 13) have demonstrated that 
the outline of the central pit crater is con- 
trolled by projectile shape and angle of in- 
cidence and that the crater depth is 
dependent on projectile density and impact 
velocity. 

Brownlee et al. (ref. 15) measured crater- 
circularities from Scanning Electro-Micro- 
scope (SEM) photographs that were taken 
with the electron-optical axis normal to the 
cratered surface. A “circularity index” was 
defined as the ratio Am/Ac, where Am is 
the area measured along the inferred inter- 
section of the surrounding target surface with 
the inside of the pit rim, while Ac is the 
area of the smallest circle which just encloses 
A,„. Circularity indices measured for 131 
micron-sized craters demonstrate the rarity 
of highly noncircular pits (fig. 3). Many of 
the noncircular craters in figure 3 are elon- 
gated and shallow, indicating that they were 
produced by oblique impact rather than 
highly irregular projectiles (refs. 15 and 16) . 
By comparison with laboratory simulations 
using irregular projectiles (ref. 14), it is con- 
cluded that highly nonspherical shapes such 
as rods or platelets are rare or nonexistent in 
the micrometeoroid complex. If dust grains 
were modeled as prolate ellipsoids, then the 
observed crater circularities suggest an aver- 
age length to width ration of < 2. 

Depth/diameter ratios were determined 
for 70 craters (ref. 15) by use of the con- 
tamination line profiling technique of Vedder 
and Lem (ref. 17) and parallax measure- 
ments from SEM stereo photos. The crater 
depth/diameter ratios refer to the maximum 
pit depth below the original uncratered sur- 
face divided by the mean diameter of the 
inside of the pit rim. Figure 4 illustrates the 
results of 70 lunar craters in histogram form 
together with laboratory cratering data of 
Vedder and Mandeville (ref. 12) . Because the 
laboratory data do not extend beyond 13 
km/s impact velocity and because the veloc- 
ity distribution of small meteoroids is not 
well known, it is not possible to determine 
exact particle densities. It is obvious, how- 
ever, that the data are entirely inconsistent 
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Figure 2 . — Typical lunar experimental microcraters on glass surfaces. Note the change of crater morphol- 
ogy with size: (a) very small lunar microcrater that displays neither concentric fractures nor a spall 
zone (note the raised, glassy rim); (h) lunar crater that displays concentric fracture zone indicative of 
incipient spallation; (c) lunar crater with partially developed spall zone; (d) lunar crater with com- 
pletely developed spall zone; (e) experimental crater (Al-projectile into soda lime glass; impact velocity: 
9.9 km/s); and (f) experimental crater (polystyrene projectile into soda lime glass; impact velocity: 5.7 
km/s. Note shallow crater depth and compare with S(a)-(e)). 




micrometeoroid complex and regolith evolution 


609 



CIUCULARITY INOiX (Am/Ac) 


Figure 3 . — Histogram of the circularity index of 131 
microcraters ranging in size from .2 to 80 ym in 
diameter (rock 15286). Though not illustrated, the 
circularity index is independent of pit crater di- 
ameter. 


with micrometeoroid densities less than unity. 
The rarity of deep craters also appears 
to exclude the possibility that a signif- 
icant fraction of particles could have den- 
sities as high as iron. Figure 4 apparently 
implies that most micrometeoroids ( < 50/jn 
in diameter) have densities in the 2 to 4 
g/cm"’ range, if one assumes an average im- 
pact velocity of 20 km/s. Even for velocities 
between 10 and 30 km/s, the above densities 
are approximately valid. 

Only « 10 percent of the total crater popu- 
lation may offer different interpretations. Of 
those exceptions, the so-called “pitless” cra- 
ters are by far the most abundant (» 80 per- 
cent). They do not possess a glass-lined pit 
(fig. 5(a)) and could be interpreted as low- 
velocity, “secondary” craters. However, as 
illustrated in figure 5(b) and as observed 
numerous times, there is strong evidence that 
many “pitless” craters did indeed have a 
glass-lined pit, which was spalled off either 
during crater formation or thereafter (ref. 
18). Thus, many of these structures are also 
potential candidates for a “primary” origin 
(ref. 19). Another exceptional crater type, 
termed “multiple pit crater,” is illustrated in 
figure 5(c) ; figure 5(d) documents a labora- 
tory equivalent produced by an agglutinate 
of minute glass spheres (ref. 12). Conse- 
quently it is conceivable that “multiple pit 
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Fig^ure 4 . — Experimentally determined depth/diame- 
ter ratios using projectiles with densities from 1 
to 7 g/cm' and impact velocities from 3 to 13 km/s. 
The inserted histogram on lunar depth/diameters 
is based on 70 craters. 

craters” are indeed caused by projectiles of 
low density and nonhomogeneous mass-dis- 
tribution, i.e., “aggregate” structure; how- 
ever, they are rare exceptions and far less 
frequent than suggested by Vemiani (ref. 
20), Hughes (ref. 21), and many others. 


CRATER POPULATIONS ON LUNAR 
ROCKS 

In analogy to large-scale lunar surfaces 
(Gault (ref. 22) ; Shoemaker et al. (ref. 1), 
and others), two basic types of crater popu- 
lations need to be distinguished: (a) “pro- 
duction” and (b) “equilibrium” populations. 
By definition, “production populations” are 
limited to rock surfaces of low, absolute 
crater densities, i.e., of short exposure pe- 
riods. With time, more and more impacts 
will occur in already cratered areas until 
finally the surface becomes so densely cra- 
tered that each new event will destroy an 
already existing one. Such a surface has 
reached “equilibrium.” “Transition popula- 
tions” are intermediate between “produc- 
tion” and “equilibrium” conditions. Most 
lunar rocks are either in “transition” or 
“equilibrium” condition; genuine “produc- 
tion populations” are rare. 

Because production surfaces exclusively 
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Figfure 5 . — Unusual craters: (a) "pitless” lunar crater (note the similarities and possible transition to crater 
5(h)); (h) similar sized lunar crater with pit (note that spallation action was severe enough to undercut 
the glass-lined pit, leaving it barely attached to the crater bottom); (c) "multiple pit" crater on lunar 
glass surface 15286); and (d) "multiple pit” crater produced in the laboratory. 


display a complete record of all craters pro- 
duced, only they are suitable to deduce mass- 
frequencies and the flux of micrometeoroids. 

Cumulative crater size distributions for 
“production populations” on samples 12054 
(ref. 4) and 60015 (ref. 23) are shown in 
figure 6; though other genuine “production 
populations” were investigated, the two 
curves illustrated are considered the best 
available over the size range indicated. The 
absolute crater densities for the two samples 
differ by almost a factor of 2, reflecting dif- 
ferent times and/or geometry of exposure. 
The relative crater size frequency, however, 
is nearly identical. 

Figure 7 illustrates “production” data 
resulting from SEM studies. The relative fre- 
quencies were normalized to surface 15205 


at a pit diameter of 1/xm. The illustrated 
data are considered the best available. The 
differences in the distributions and the pres- 
ence of an inflection at pit diameters between 
ljum and 10/tm are subject to a variety of 
interpretations. They will be discussed later. 

Because the rock surfaces that have 
reached “transition” and/or “equilibrium” 
conditions are less suitable for study of the 
micrometeoroid complex, they will not be 
treated extensively here (refs. 5, 23, 25, 26, 
and 27). However — if coupled with solar 
flare track exposure ages — ^they may still 
contribute to the flux determination of mi- 
crometeoroids ; minimum fluxes may be ob- 
tained, because a number of the craters 
produced are destroyed and not observable 
anymore. 
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PIT DIAMETER, MICRONS 

Figure 6 . — Typical binocular crater-size frequency 
distributions for lunar glass surfaces in production 
state (1205i is based on 960 craters; 60015 is 
based on 665 craters). 



Fig^ure. 7 . — Typical scanning electron microprobe 
crater-size frequency distributions for small micro- 
craters on lunar glass surfaces in production state 
(15205 is based on 1100 craters, 15286 on ^ 500, 
and 15017 on •=> 300). 


MASS-FREQUENCY OF 
MICROMETEOROIDS 

Crater simulation experiments provide the 
only basis to obtain information concerning 
the mass distribution of micrometeoroids by 
converting crater dimensions into projectile 
parameters. The physical processes governing 
impact cratering are complex and presently 
not understood in great detail, despite con- 
siderable laboratory work. Especially, the 
energy partitioning for small- and large-scale 
cratering and the effects of target strength, 
gravitational forces, and varying impact 
velocities, i.e., appropriate “scaling laws,” are 
still subject to experimental work that ulti- 
mately will result in a theoretical under- 
standing. Therefore, extrapolations from 
laboratory data may allow the use of a 


variety of empirical calibration approaches. 

Four basic calibration techniques for 
microcraters are currently in use (fig. 8). 
Two are based on electrostatic dust acceler- 
ator experiments (refs. 6, 9, 24, and 30), and 
two calibration techniques utilize results 
from ballistic ranges (refs. 28 and 31), 
while Nagel (ref. 29) employed a lithium 
plasma gun (for more detail, see ref. 16). 

Relative crater-size frequency distributions 
ranging from .1 to almost 1000 microns in pit 
diameter may be constructed from the data 
presented in figures 6 and 7 by normalizing 
the absolute crater densities with respect to 
exposure time, exposure geometry, and sur- 
face area. An important assumption under- 
lying such a normalization is that these 
relative frequencies remained constant with 
time, because surfaces of different crater 
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CALIBRATION CURVES 



Figfure 8 . — Various calibration methods presently in 
use to derive micrometeoroid masses from mea- 
sured pit diameters (Dp) or spall diameters (D.). 
The ratio D,/Dp is variable from rock to rock with 
values between S.8 and U-5 on lunar glasses. Note 
that agreement between various techniques is close, 
if a D,/Dp of i.S is applied. 

Log ,0 I MPACT ENERGY, ergs 



Figure 9 . — Differential frequency of pit diameters 
and their corresponding particle mass and energy 
distributions. The binocular data (ISOSi, 60015) 
and SEM data (15205) are joined at a pit diameter 
of 100 pm as indicated. 


densities, i.e., different absolute exposure 
times, need to be normalized. Figure 9 shows 
such a normalized, differential crater fre- 
quency distribution based on glass-surfaces 
12054, 60015, and 15205. The corresponding 
mass- and energy-scales are based on the cal- 
ibration by Gault (ref. 28) as shown in 
figure 8. For masses > 10"*®g ( = impact en- 
ergies above 200 ergs) this distribution is 
in basic agreement with that obtained by 
satellite- and ground-based measurements 
(refs. 32 and 33). Though the irregularity of 
the distribution at lower masses will be more 
thoroughly discussed later, it can already be 
seen that 

1. Particles in the 10“*® to 10**® g range 
are most numerous. 

2. The bulk of the meteoroid mass or en- 
ergy impacting the Moon is confined to par- 
ticles 10*® to 10*^‘ g in mass (see also refs. 
4 and 34). 

FLUX OF MICROMETEOROIDS 

Micrometeoroid fluxes are obtained by cor- 
relating absolute crater densities with the 
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Figure 10 . — Correlation of surface residence times 
of specific rocks mostly determined with solar flare 
tracks (see table 1) and absolute frequency of pit 
craters > 500 pm in diameter/cm'. A straight line 
going through the origin establishes the crater pro- 
duction rate, i.e., the flux of micrometeoroids 
> 2 X JO"' g. A best estimate for the average flux 
over the past 10" years is ^ 5 pits/cm'. 





MICROMETEOROID COMPLEX AND REGOLITH EVOLUTION 


613 


absolute exposure age. A summary of such 
correlations for binocular crater counts on 
selected rocks is illustrated in figure 10, with 
the cumulative crater frequency for pits 
above SOO/^m in diameter. Most data points 
shown lie below possible correlation lines 
and therefore are in or approaching equilib- 
rium with respect to cratering. A correlation 
line corresponding to a crater production 
rate of five pits with diameters equal to or 
greater than 500 microns per cm^ per mil- 
lion years lies within a factor of 2 of data 
for 12054, 12017, 12038, and 14301. Upon 
visual inspection of these samples, only rock 
12038 was not clearly in production with 
respect to cratering. A factor of 2 is the 
estimated uncertainity in the solar flare 
track method used for the exposure time 
measurements. 

Another approach to measure the meteor- 
oid flux and possible changes with time has 
been pursued by Hartung et al. (ref. 35). 
Separate solar flare track exposure ages were 
determined for 56 individual pit craters 
larger than 20/tm on rock 15205. The results 
illustrated in figure 11 indicate that the for- 
mation ages of these craters are not uni- 


formly distributed; significantly more cra- 
ters are produced during the last 10 000 
years. Thus it appears that the present-day 
micrometeoroid flux is enhanced over that of 
the past 10^ to 10® years by slightly more 
than an order of magnitude. The values ob- 
tained for the past 3000 years are in good 
agreement with present-day satellite mea- 
surements (refs. 33 and 34). 


Discussion of the Micrometeoroid 
Complex 

IMPLICATIONS OF PHYSICAL 
PROPERTIES 

On the basis of laboratory cratering ex- 
periments, the morphologies of microcraters 
are interpreted to indicate that they were 
formed by equidimensional, nonporous pro- 
jectiles of densities between 2 and 4 g/cm®, 
which impacted with velocities in excess of 
5 km/s. These results are in part contrary 
to popular hypotheses and they may have 
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Fi^re 11 . — Preliminary formation ages of 56 individual microeraters ranging in size from 20 to SOO mi- 
crons in pit diameter on glass surface 15205,51. Note the steep increase in crater production rate between 
0 and 10' years, corresponding to a twofold increase every 3000 years. 
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some significant astronomical consequences. 

A cometary origin for micrometeoroids is 
strongly suggested by a variety of indepen- 
dent analyses (e.g., refs. 33 and 36). The 
particulate matter within comets is believed 
to represent unfractionated, solar abundances 
similar to Type 1 carbonaceous (Cl) chon- 
drites. Our mass densities are entirely con- 
sistent with Cl chondrites, the constituents 
of which range in density from approximately 
1.5 g/cm"* for aggregates of phyllo-silicates 
to magnetite grains of Sg/cm’ in density 
(ref. 37). 

Much lower densities with an average of 
.5 to .8 g/cm"* have been suggested for the 
somewhat larger meteors, i.e., particles 
> 1(H g (refs. 20 and 21). Though our de- 
tailed analysis of crater morphologies is con- 
fined to craters below lOO/^m in diameter, 
i.e., particles < 10“* g, even pit craters larger 
than 1 cm, caused by particles of approxi- 
mately 10-® g, display qualitatively the same 
morphologies. Though precise laboratory cal- 
ibrations are not available for such large 
structures, we suggest that most particles of 
10-® to 10~® g may also have a density of 
more than unity. 

The equidimensional character of micro- 
meteoroids may also have significant astro- 
nomical implications, if we accept a cometary 
source. Traditionally it is suggested that such 
materials are similar if not identical in 
chemistry and shape to grains found in car- 
bonaceous chondrites, because they are be- 
lieved to represent primordial condensates 
from similar environments in the solar neb- 
ula. These grains are thought to be vapor 
growth products of highly nonspherical 
shape like platelets, rods, and whiskers (refs. 
14, 38, 39, and 40). Such grains were ob- 
served in a variety of carbonaceous meteor- 
ites, e.g., Allende, which is thought to be a 
fine example of “early condensates” ( ref. 41 ) . 
Clearly the microcrater circularities are in- 
compatible with such elongated grains. These 
findings either imply that the postulated 
grain shapes are incorrect and virtually non- 
existing in the environment of comet forma- 
tion or that the micrometeoroid complex is 
also the result of multiple collisional events 


prior (!) to incorporation into cometary 
matrices. Recent developments in meteorite 
research provided strong evidence that col- 
lisional processes in the early history of the 
solar system may have played a dominant 
role. 

Regardless of what caused the microme- 
teoroids’ equidimensional if not spherical 
shape, needles, platelets, rods, whiskers, and 
other elongated or irregular particles seem 
not to make up a significant part of cometary 
silicates, if one accepts at all a cometary 
source area. The possibility that most of 
these particles constitute debris of collisional 
processes during accretion, rather than pri- 
mary condensates, cannot be excluded. 

MASS-FREQUENCY 

The frequencies of micrometeoroid masses 
ranging from 10"*® to 10-'* g are summarized 
in figure 12, together with a variety of satel- 
lite- and earth-based measurements. Two 
types of microcrater frequencies are ob- 
served: that displayed by samples 15205, 
15076, and 15017 ; and that of sample 1^286. 
Though experimental conditions (most dom- 
inantly target-smoothness and total number 
of craters counted) may be responsible for 
subtle differences of the first type, the dif- 
ferent behavior of 15286 seems beyond sta- 
tistical error. Rock 15205 is based on 950 
craters and sample 15286 on 500 craters. 
Thus, two questions remain: (1) Why are 
there two different frequency types? and (2) 
What causes an apparent bimodal mass dis- 
tribution? 

Sample 15286 is unique, though there are 
other samples (e.g., 12024,81 and 14257,F 
(ref. 30)) that may be similar. Their dif- 
ferent mass-frequencies may be caused by 
extreme solid angles of exposure (ref. 23) 
that effectively influence the energy distri- 
bution, because of the increased effects of 
oblique impact (ref. 28). It is also conceivable 
that such surfaces were essentially pointing 
toward lunar North, i.e., out of the ecliptic 
plane, where they potentially could intercept 
a different population of cosmic dust than 
within the ecliptic plane. 


MICROMETEOROID COMPLEX AND REGOLITH EVOLUTION 


615 



Figure 12 . — Comparison of the cumulative mass-frequency slopes of a variety of observational techniques, 
but, in particular, of five 'well-documented lunar glass surfaces. Most individual satellite data do not give 
a differential flux; the position of the satellite data was constructed by extrapolating the slope from the 
cummulative mass-frequency curve of Dohnanyi (ref, S3, fig. 1). The length of the bars indicates the mass- 
range over which the corresponding slope is valid. 


Curves 15205, 15076, and 15017 are be- 
lieved to be typical for micrometeoroids im- 
pacting the Moon, simply because such 
distributions are the most frequent ones. 
Samples 60502,17; 15927,3; 15301,79 (ref. 
24) ; and 15015 (ref. 5) yield similar results. 
The cause of this apparent bimodal mass dis- 
tribution is presently unknown. However, it 
is conceivable that the larger masses repre- 
sent the cometary particle population that is 
spiraling totvard the Sun. During and upon 
solar approach, individual particles may suf- 
fer fragmentation as well as melting and/or 
vaporization; both processes would result in 
numerous particles of very small sizes. Upon 
close solar approach they may be propelled 
atvay from the Sun again by solar radiation 
and have a second opportunity to encounter 
the lunar surface (refs. 36 and 42). 

Subtle differences in crater populations 
may yet be caused by a completely different 


mechanism. Morrison et al. (ref. 5) and 
Blanford et al. (ref. 43) report that lunar 
rock surfaces are significantly modified on the 
micron scale by the accretion of regolith 
particles; most dominantly, disk-shaped, 
glassy splashes and droplets. These accre- 
tionary objects are so numerous that they 
accumulate obviously at a faster rate than 
the surface is destroyed by microcraters. 
Given sufficient time they may even build up 
layers of a few microns in thickness, giving 
some of the hand specimens a typical, pati- 
nated appearance. Thus a “constructive” 
accretion process is competing with the 
“destructive” cratering process and the 
micron-size crater population may be some- 
what modified. The unambiguous presence of 
particles below 10"^® g in mass, however, 
negates the existence of a radiation pressure 
cutoff. According to Gindilis et al. (ref. 44), 
the lack of such a cutoff is highly compatible 


616 


COSMOCHEMISTRY OF THE MOON AND PLANETS 


with particle densities of 2 to 4 g/cm®, i.e., 
with silicates, for which gravitational forces 
appear to dominate radiation pressure. This 
result corroborates our conclusions about 
particle densities. 

Figure 12 also illustrates one fundamental 
advantage of lunar glass-surfaces as micro- 
meteoroid detectors: at present, the lunar 
rock detector spans 12 orders of magnitude 
in mass and thus possesses a “dynamic range” 
duplicated nowhere. The potential identifica- 
tion of a bimodal size distribution is only due 
to such a large dynamic range. 

Additional work with carefully selected 
samples is required to clarify what causes the 
two basic frequency types and the apparent 
bimodal distributions. The above explanations 
have to remain tentative until carefully selec- 
ted and precisely oriented surfaces are inves- 
tigated in detail. 

MICROMETEOROID FLUX 

A detailed comparison of micrometeoroid 
fluxes derived from lunar sample analyses 
and satellite measurements is presented in 
figure 13. It is impossible to discuss each de- 
tail and thus we offer a few general com- 
ments only, quoting Horz et al. (ref. 16) : 

The moon is a rotating sampler, and 
the directional distribution of micro- 
meteoroids is extremely non-uniform as 
shown by Berg and Griin (1973) (ref. 
48) and Hoffmann et al. (1973) (ref. 
47). Accordingly, the meteoroid flux 
differs about 3 orders of magnitude be- 
tween the direction of the earth’s apex 
and anti-apex. Furthermore, particles 
> 10-^2 g are collected almost exclusively 
during the apex orientation of the Pio- 
neer and HEOS sensors. Hence, in this 
mass range, also the Moon may collect 
particles from only the apex direction. 

As a consequence, a “detector” on the 
rotating lunar surface can “register” 
meteoroid impacts effectively only part 
of the time. Therefore, fluxes derived 
from lunar crater statistics may have to 
be increased by as much as a factor of 
7T for comparison with satellite data that 


were taken in the apex direction. Also, 
apex-pointing data generally have been 
corrected upward to a standard 27r-sterad 
exposure angle, assuming an isotropic 
flux. Thus, an actual anisotropy (as re- 
ported by the HEOS and Pioneer ex- 
periments) leads to an overestimation of 
the flux. Therefore, the satellite results 
seem to represent an upper limit for the 
flux. 

The “apex” particles show an average 
impact velocity of only 8 km/sec (Hoff- 
mann et al., 1973) (ref. 47). The fluxes 
from lunar rocks, however, are calculated 
with a standard velocity of 20 km/sec. 
The necessary corrections will increase 
the projectile masses and thereby effec- 
tively enhance the Moon-based flux for 
masses > 10“*" g by a factor of approxi- 
mately 5. 

The situation for masses < 10“** g is 
highly complex. Berg and Griin (1973) 
(ref. 48) have reported that most events 
of these masses occur with particles that 
have relative velocities of at least 50 
km/sec. The lunar flux curves given for 
these masses in fig. 12 are, however, 
based on a 20 km/sec impact velocity; if 
corrected to 50 km/sec, they will shift 
towards smaller masses, possibly as 
much as a factor of 10. 

As a consequence, the fluxes derived from 
lunar crater statistics may agree within the 
order of magnitude with direct satellite re- 
sults if the above uncertainties in velocity 
and directional distribution are considered. 

Figure 14 presents some basic constraints 
derived from a variety of independent lunar 
studies on the flux of micrometeoroids and 
larger objects. The only direct measurements 
are the impact events registered by the Pas- 
sive Seismic Experiment (ref. 62) and the 
micrometeoroids encountered by the space- 
craft windows (ref. 57). Upper limits on the 
flux can be derived from more cratering rate 
(refs. 63, 64, and 65). Accordingly, the 
flux over the past 3.0 X 10® years has re- 
mained fairly constant. The “geochemical” 
evidence is based on the abundance of sider- 
ophile trace elements indicative of type and 
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O (ref. 49) 
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▼ (ref. 51) 
V (ref. 52) 
<0 (ref. 52) 

<C> (ref. 53) 
^ (ref. 54) 
♦ (ref. 55) 
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Moon (Hartung et al., 1972b) (ref. 57) 

Moon (Morrison et al., 1973) (ref. 58) 

Moon (Schneider et al., 1973) (ref. 59) 

Moon (Neukum et al., 1972) (ref. 59) 

OGO II, IV (Nilsson et al., 1969) (ref. 60) 

OGO III (Alexander et al., 1972) ^ef. 61) 


HEOS 2 (Hoffmann et al., 1973, apex) 
Pioneer 8 (Berg and Grun, 1973) 

Mariner II (Alexander et al., 1965) 

Mariner IV (Alexander et al., 1965) 

Various Satellites (Nazarova, 1968) 

Lunar Explorer 35 (Alexander et al., 1971) 
Explorer XXIll (Naumann et al., 1969) 
Pegasus I, II, III (Naumann et al., 1969) 
Explorer XXIll (O’Neal, 1968) 

Cosmos 135 (Konstantinov et al., 1969) 
Cosmos 163 (Konstantinov et al., 1968) 
Lunar Orbiter 1-5 (Gurtler and Grew, 1971) 


NASA Model (Cour-Palais, 1974) 
Hawkins, 1963 (photographic) 

Moon, prehistoric flux (Neukum, 1973) 
Moon, historic flux (Neukum, 1973) 
(Radar) Kaiser, 1961 
(Visual) Lindblad, 1967 


Figure 13 . — Comparison of lunar and satellite micrometeoroid flux data. 
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Figure 14 . — Constraints on the flux of micromete- 
oroids and larger objects, according to a variety of 
independent lunar studies. 


amount of meteoritic contamination in the 
lunar soil (ref. 66). Erosion rates on lunar 
rocks range from approximately .2 to 2 mm/ 
10* years (refs. 67, 68, and 69). Taking the 
highest erosion rate and applying cratering 
data of Gault (ref. 28) make possible the 
definition of an upper flux limit. Further- 
more, the negative findings on the Surveyor 
III camera lens (ref. 70) and the perfect 
preservation of the footpad print of Surveyor 
III (ref. 71) also define an upper limit. A 
lower limit results from the study of solar 
and galactic radiation tracks in lunar soils 
(refs. 72, 73, and 74). It is found that some 
cm-thick layers of regolith have resided on 
the lunar surface essentially undisturbed for 
1 to 2 X 10^ years. Because the regolith 
is believed to be reworked by micrometeor- 
oids only, the flux could not have been sig- 
nificantly lower than indicated; otherwise 
still older residence times for the soil layers 
would be obtained. Strictly, only the passive 
seismometer, the Apollo windows, and the 
mare craters yield a cumulative mass distri- 
bution. All other parameters are a bulk mea- 
sure of either meteoroid mass or energy ; the 
corresponding “flux” was calculated using 
the differential mass distribution obtained 
from lunar microcraters. Accordingly, the 
corresponding arrows may be shifted any- 


where along the line defining the “upper” and 
“lower” limits. 

The data shown in figure 11 (ref. 75) sug- 
gest that the present flux is significantly 
higher than the average flux over the last 
10* to 10® years (ref. 35). Gault et al. (ref. 
34) and Morrison et al. (ref. 5) were the 
first to indicate such a possibility because 
absolute lunar rock exposure ages, erosion 
rates, and survival times of rocks appeared 
to be incompatible with computed values that 
were based on present-day meteorite fluxes 
derived from satellites. Neukum (ref. 59) ex- 
panded on these interpretations and his 
“historic” and “prehistoric” fluxes are incor- 
porated in figure 13. Because the annealing 
behavior for radiation tracks during long- 
term exposure in the lunar environment is 
not well known and because all potential er- 
rors — both in the age dating as well as 
crater-counting — enter these considerations, 
a “historic” and “prehistoric” flux can only 
be tentatively proposed at present. The data 
of Hartung et al. (ref. 35) present the strong- 
est evidence to date. 

Although the magnitude of the flux may 
have varied over geological times, the mass- 
frequency distribution appears to have re- 
mained fairly constant. Frequencies mea- 
sured on surfaces that were constituents of 
the soil (15927, 15301, 15001, and 60502) 
most likely do reflect the meteoroid bombard- 
ment that is older than that of most rocks. 
Their size-frequency distributions agree 
within the accuracy of measurement with 
“recent” crater populations. 

Brownlee and Rajan (ref. 76) and Rajan 
et al. (ref. 77) discovered microcraters that 
are identical to lunar craters on the surface 
of glassy spherules, dislodged from the 
interior of the Kapoeta meteorite. This mete- 
orite is a loosely consolidated microbreccia 
and a striking meteoritic analog to lunar soil 
breccias in many aspects. The formation age 
of Kapoeta is approximately 4 X 10® years 
(Rajan, personal communication, 1974). 
Within the counting accuracy, the size-fre- 
quency distribution of the Kapoeta micro- 
craters is identical to lunar ones. Brownlee 
and Rajan furthermore dated one spherule 
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via solar flare tracks and derived a micro- 
meteoroid flux that is within an order of 
magnitude of the present-day flux. Because 
the track retention over 4 X 10® years in 
glassy materials is poorly known, however, 
this exposure age and the resulting micro- 
meteoroid flux has still larger uncertainties 
than young lunar glass surfaces. Blanford et 
al. (ref. 43) report numerous microcraters on 
feldspars separated from the very bottom 
(soil-sample 15001) of the 240-cm-long 
Apollo 15 drill core. The observed crater size- 
frequency distributions are essentially iden- 
tical to those of rock samples 15205 and 15017 
(see fig. 12). Because this soil was deposited 
at its site of collection more than 400 m.y. 
ago (ref. 78) the observed crater populations 
must have formed prior to that time. Micro- 
meteoroid craters are also found in virtually 
every “soil-breccia” as well as genuine soil 
samples, though their actual geological time 
period of exposure is not known at present 
(ref. 24). Taking typical noble-gas exposure 
ages of lunar soils as statistically representa- 
tive average values of the individual compo- 
nents, it may safely be concluded that 
micrometeoroid bombardment was active 
throughout geological time. From the pres- 
ently available microcrater size-frequency 
distributions it also may be concluded that 
the mass-frequency distribution of microme- 
teoroids has not changed significantly if at 
all. 

The studies on surfaces of old exposure 
ages demonstrate another unique character- 
istic of the “lunar rock micrometeoroid 
detector”: it is principally possible to de- 
lineate the flux and potential variations 
thereof through geologic history. Such poten- 
tial variations are of considerable interest to 
an explanation of the formation of the solar 
system for a variety of reasons: 

1. The presence of a minimum micromete- 
oroid mass may be determined as a 
function of geologic time. This mass, in 
turn, may be used to calculate upper 
limits on the solar radiation pressure 
and thus to the luminosity of the Sun. 
Brownlee and Rajan (ref. 76) have at- 
tempted such calculations based on the 


minimum crater diameter observed on 
the Kapoeta materials and they con- 
cluded that the solar luminosity at 
<=» 4 X 10* years was not higher than 
1.7 times its present value. 

2. The main source of micrometeoroids 
has to be sought in short-period comets. 
Significant variations in the flux of me- 
teoroids may be related to short-period 
comet “activities,” i.e., to an uneven, 
possibly sporadic rate of comet en- 
counters that are capable of putting 
micrometeoroids with bound orbits into 
the inner solar system. In addition to 
these relatively short-term fluctuations 
(millions of years) it is also possible 
that the rate of comet injection into the 
inner solar system has undergone a 
long-term secular change due to a 
general depletion of the comet inven- 
tory. 

3. Micrometeoroid detectors onboard Pio- 
neer 8/9 have intercepted a non-neg- 
ligible fraction of interplanetary par- 
ticles that have hyperbolic orbits and 
thus are interpreted to be of inter- 
stellar origin (ref. 48). Thus lunar 
rocks offer a potential opportunity to 
study interstellar grains. 

Most of the above possibilities, however, 
will require substantial amounts of work and 
are — at present — considered exciting chal- 
lenges for future research. They are men- 
tioned above only to stress the uniqueness and 
exciting potential of cratered lunar rock sur- 
faces. 

Lunar Regolith Dynamics 

The lunar regolith is a layer of fragmental 
debris of variable thickness that lies upon 
fractured bedrock. Photogeologic investiga- 
tions and detailed analysis of returned lunar 
materials revealed that repetitive meteoroid 
bombardment has been responsible for the 
formation of this layer to such an extent that 
other geological processes may be excluded. 
Impact cratering controls the overall growth 
of regolith; the lateral and vertical redistri- 
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bution of material ; the downslope mass 
wasting ; the mixing and degree of homogeni- 
zation of individual layers; the erosion of 
lunar rocks; the evolution of regolith grain 
sizes; the formation of impact melts, agglu- 
tinates, and breccias'; the migration of vola- 
tile elements ; and the admixture of meteoritic 
components and other parameters that make 
up the physical, chemical, and petrographic 
characteristics of lunar “soils.” As a conse- 
quence, it appears appropriate to combine 
observational lunar crater data and experi- 
mental impact crater mechanics into compu- 
tational models to arrive at a theoretical 
understanding of these processes. 

A variety of computational results concern- 
ing mass-movement, erosion rate of rocks, 
etc., are available (refs. 34, 59, 63, 79, and 
others). However, all these analyses suffer 
from the fact that they yielded only “average” 
values because the computations did not ac- 
count for the vagarities of the random impact 
process. Models that do, however, account for 
the randomness of the impact process both 
in space and time have been developed re- 
cently and are described below. The models 
may be used to gain a qualitative if not quan- 
titative insight into some of the above pro- 
cesses. Some of these models consider craters 
up to 1500 m in diameter and thus are of 
drastically different dimensions than the 
craters treated in the preceding sections. 
Furthermore, it is also important to note that 
the models are principally independent of the 
absolute flux of meteoroids. The time param- 
eter is linearly related to the total number 
of craters produced. Thus model elapsed times 
can easily be converted into absolute times 
by applying the best estimate of the absolute 
meteoroid and micrometeoroid infall rates. 


LARGE-SCALE REGOLITH CRATERING 

The gross-accumulation of the regolith 
debris layer has been the subject of a variety 
of treatments, e.g., Marcus (ref. 80) and 
Shoemaker (ref. 63). It has been demon- 
strated that the overall regolith thickness 
increases with increasing numbers of craters 


that range roughly in diameter from 10 to 
1000 m. Oberbeck and Quaide (ref. 81) 
pointed out that the growing debris layer acts 
as a buffering medium and thus strongly con- 
trols the geometry of different crater sizes. 
Accordingly, the actual thickness for a given 
lunar surface area can be related to the total 
number of craters produced as well as to the 
relative frequencies of differently shaped 
craters such as “normal,” “flat-bottomed,” 
“concentric,” and “central mound” craters. 

Oberbeck et al. (ref. 82) have developed 
a large-scale Monte Carlo-based computer 
program that simulates the evolution of the 
regolith and that also predicts the relative 
frequencies of the above four basic crater 
morphologies for any given regolith thick- 
ness. It is important to note that these cal- 
culations were performed with observed 
lunar cratering parameters, i.e., detailed 
crater geometries and distributions of asso- 
ciated ejecta blankets. No cratering scaling 
laws needed to be assumed. 

A crater production size-frequency distri- 
bution of N = KD-^ * was empirically deter- 
mined and used throughout these calculations 
{N = cumulative number of craters larger 
than diameter D, i.e., > 1 m) . Some pertinent 
results are discussed below; for detailed in- 
formation the reader is referred to Oberbeck 
et al. (ref. 82). 

Figure 15 illustrates the relationship of 
the calculated median regolith thickness 
(Rm) as a function of absolute numbers of 
craters produced. A relationship of 

R„ = 6.2 X 10-®A« (1) 

is derived and may be used to predict the 
median thickness for any surface area where 
crater size-frequency distributions can be de- 
termined and where the cumulative crater 
production distribution has the form of N — 
K • 

However, the regolith thickness is variable 
over distances measured in hundreds of me- 
ters as evidenced by high-resolution photog- 
raphy and field inspection by the astronauts, 
despite the fact that the overall reference 
surface must have been exposed to the mete- 
oroid bombardment for the same period of 
time. Figure 16 compares actually measured 
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Figure 15 . — The overall regolith growth as a func- 
tion of craters produced, i.e., time ( see equation 1 ). 


thickness distributions (ref. 81) with those 
obtained in the Monte Carlo simulations. 
The agreement is good and lends additional 
support to the hypothesis that the regolith at 
the sites investigated by Quaide and Oberbeck 
(ref. 81) is primarily caused by impact com- 
minution processes. 

However, the above Monte Carlo model on 
regolith formation yielded additional infor- 
mation: with increasing thickness of regolith, 
only larger and larger craters are capable of 
penetrating the existing, buffering debris 
layer. Thus, with increasing time, it takes 
larger and larger craters to excavate pristine 
bedrock. The Monte Carlo simulations there- 
fore continuously monitored, per each crater 
size class, the total volume excavated from the 
pristine substrate (F,) and the already ex- 
isting regolith layer (Vr) throughout the 
time required to build up the regolith to a 
given thickness. Figure 17 illustrates the ra- 
tio V,/Vr for three different regolith 
depths. The ratio YJVr is a function of 
crater diameter and is described by 

F./F« = C‘ D" (2) 

where C is a constant for a given distribution 
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Figure 16 . — Empirically determined regolith thick- 
ness distribution for four different lunar surfaces 
measuring >=> 200 km’ each. The empirical determi- 
nation is based on the abundance of various crater 
geometries reflecting the presence of a competent 
substrate. Note the good agreement between ob- 
servations and Monte Carlo cratering simulations. 


of craters (n = 1 — 1.3) ; furthermore, C can 
be related to K in the crater distribution 
expression N = KD~^ * by 

C = 1.02 X (3) 

and by substitution 

V./Vr = 1.02 X (4) 

Thus, over the range of values of K charac- 
teristic for mare terrains (as an example), 
the effective size boundary between mixing 
and new debris-producing craters becomes 
progressively larger. The average mixing 
zone therefore becomes deeper. Accordingly, 
older and thicker reolith deposits should be 
more thoroughly reworked than more youth- 
ful ones. 

Figure 18 illustrates the cumulative con- 
tributions of various-sized craters that have 
built up a regolith layer of 4.7 m in median 
thickness. It is obvious from figures 17 and 
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Figure 17 . — Relative contributions of pristine bed- 
rock from the “substrate” for various crater sizes 
and regolith depths. The “volume regolith” is that 
volume that is reworked debris excavated by prior 
cratering. Note that predominantly the larger 
craters excavate bedrock and thus chiefly contrib- 
ute to the overall regolith growth with increasing 
regolith thickness, i.e., time. 


18 that relatively small craters (e.g., < 10 m 
in diameter) have contributed significant 
amounts to the overall regolith, but it is also 
readily seen that these contributions occurred 
while the regolith was relatively thin, i.e., in 
the early stages of regolith formation. At 
present it is predominantly structures > 100 
m in diameter that control the overall rego- 
lith growth while the smaller structures are 
confined to reworking these materials. As 
a consequence, the regolith thickness in- 
creases in general and in particular during 
its more recent history (i.e., the past « 10® 
years) due to the effects of relatively large 
cratering events that are capable of exca- 
■vating pristine bedrock. This newly added 
material will always be delivered on top of 
the existing debris in discrete swaths of 
ejecta. The regolith therefore has to be envi- 
sioned as a complex sequence of numerous. 



Figure 18 . — Contributions ( = volume percent) of 
various source areas at depth “d” to the overall 
composition of a typical mare regolith having a 
median thickness of U.7 m. 


overlapping ejecta blankets. These discreet 
blankets constitute some first-order disconti- 
nuities and heterogeneities in the evolving 
regolith. We will demonstrate in the next 
chapter that it is principally possible to pre- 
serve parts of these blankets despite heavy 
meteorite bombardment. Though there will 
be extensive mixing, there will not be com- 
plete homogenization of the regolith. 

SMALL-SCALE REGOLITH CRATERING 

It is obvious from the preserved stratig- 
raphy in returned core tube samples that 
reworking has not obliterated all stratifica- 
tion in the regolith. It is just as obvious, 
however, that every stratum that resided at 
the very lunar surface has been subjected to 
the meteoroid bombardment and the rework- 
ing process which — due to the mass- 
frequency distribution of interplanetary 
matter — does operate on a micron-to-meter 
scale. The extent to which a stratum survives 
thus must be a function of its original thick- 
ness and length of surface residence time be- 
fore it is blanketed by ejecta of sufficient 
thickness to effectively remove it from the 
active zone of reworking. Absolute parame- 
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ters for these variables principally vary with 
absolute time ; i.e., the cumulative number of 
craters produced. The absolute number of 
craters that contributed to the history of re- 
turned samples must certainly be larger than 
the numbers presently observable in the re- 
spective sampling areas because these are in 
crater-saturation for craters < 100m in 
diameter (refs. 22 and 63). Thus the poten- 
tial surface history of sampled materials can 
be understood only if a continuous bombard- 
ment history is assumed in computational 
models. 

Because meteoritic impact is a random 
process, any given point on the lunar surface 
has a unique history as compared with any 
other given point. On the other hand, the 
dominant role of meteoroid impact suggests 
that over extended periods of time any two 
areas of a given size will have experienced 
similar histories that differ only in details to 
a greater or lesser degree. Thus computa- 
tional analyses that yield “average” values 
may be useful in understanding the basic 
processes. However, they should be applied 
only with extreme caution to actual sample 
data because of the uniqueness of each indi- 
vidual sampling location. “Averages” are cer- 


tainly a valid framework for returned sample 
interpretations; but they should be applied 
only if sufficient statistical sample data are 
available. For individual data point such 
averages cannot be applied and may lead to 
grossly erroneous results, because significant 
deviations from the “average” have to be ex- 
pected from a random process. 

Gault et al. (ref. 83) have shown that the 
probability P„ of a given point on the lunar 
surface remaining undisturbed, i.e., lying 
outside a crater of apparent diameter D in a 
time interval t is given by 

= exp. ( - xATtZ)V4) (5) 

where N is the flux of the randomly distrib- 
uted impacting bodies per unit time and 
area which produce craters of diameter D. 
The probability Pc of a given point’s having 
been affected, i.e., lying within exactly n cra- 
ters of size D can be expressed as 

Pc(n)=Pu{^NtDyAr/,{\) ( 6 ) 

Equation 6 is the Poisson probability func- 
tion. Using the values given by Molina (ref. 
84) for a range of n = 0 — 153 and 
(ttNID^/A) = .001 — 100 and calculating 
additional terms up to n = 10®, Gault et al. 
(ref. 83) calculated how many times a given 
surface area may be impacted. A microme- 



Figfure 19 . — Analytical model baaed on Poisson probability function describing the amoynt of surface area 
(percent) that will be affected by meteoriod impact, and the number of times impact will have occurred 
after given model elapsed times. 
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teoroid mass-distribution of the form N = 
1.45 was used for 10"“ to lO-'^g me- 
teoroid mass (to) and N — 9.14 X 10 ®to* *i® 
for projectiles 10-^ to 10® g. Furthermore 
a standard impact velocity of 20 km/s, 
together with laboratory cratering data into 
unconsolidated materials (ref. 28), was ap- 
plied in these calculations. The principal re- 
sult is shown in figure 19. 

Virtually identical results (fig. 20) were 
obtained in a Monte Carlo-based computer- 
simulation by Horz et al. (ref. 85), that ap- 
plied the crater size-frequencies of figure 6 
and a random number generator to determine 
impact coordinates and the magnitude of 

PROBABILITY (%) 
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Figrure 20 . — Almost identical data as in figure 19, 
resulting from a Monte Carlo computer program. 
Note the multiple bombardment history of frac- 
tional surface areas ivith increasing time, Le., 10‘ 
craters produced. 


each cratering event. The curve labeled “lx” 
in figure 20 indicates how much surface area 
is affected at least one time. Note that 50 per- 
cent of the test surface ( = 44 cm®) is already 
cratered after 8300 craters, 152 to 22500 ^m 
in spall diameter. It takes more than a factor 
of 10 additional craters to affect the remain- 
ing 50 percent. Though qualitatively not 
surprising, these absolute numbers were 
unexpected. Furthermore, figure 20, for ex- 
ample, illustrates that by the time 99 percent 
of the surface is cratered at least IX (99 
percent probability), 92 percent of the sur- 
face is already cratered twice, 81 percent has 
suffered at least three impacts, 59 percent is 
cratered four times, etc. As 99.99 percent of 
the surface is cratered at least once, 88 per- 
cent will already be affected at least five 
times, etc. 

An extension of the data illustrated in 
figure 20 is presented in figure 21, which is 
based on 10® craters (ref. 85). Per each 
model-elapsed time, it was determined how 
often a given fractional surface area was im- 
pacted. Note that when the entire area 
( = 100 percent) is cratered at least one time, 
50 percent has suffered already 12 impacts 
and 10 percent surface was cratered at least 
17 times. Or, alternatively, if it takes one 
time to affect 50 percent of a lunar surface, 
it will take 3.8 times longer to affect 90 per- 
cent, a factor of 6.6 longer to cover 99 per- 



' 10 100 
MODEL TIME 


Figure 21 . — General probability of multiple bom- 
bardment history for various fractional surface 
areas (total craters produced: 10' ; see text). 
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cent and, finally, 19 times longer to crater 
100 percent of the surface. The model times 
indicated in figures 19 through 21 will be 
used in identical fashion throughout this 
report. Unit time is defined as the time re- 
quired to affect 50 percent of the surface 
area at least one time. 

Figures 19 through 21 illustrate a funda- 
mental characteristic of the impact process. 
While finite — though admittedly small — 
surface areas may remain unaffected for 
long time periods, other areas have already 
suffered repetitive bombardment. Conse- 
quently, within any cratered terrain, small 
surface areas may be encountered that have 
dramatically different bombardment histo- 
ries despite the fact that they were exposed 
to the same micrometeoroid environment for 
the same period of time. 

We now turn to the mixing of the regolith. 
The above models are also a measure of how 
much kinetic energy is deposited randomly 
in space and time into a unit area of lunar 
surface. Therefore, one can associate with 
that energy either a crater diameter (as 
above) or a corresponding crater depth. 
Gault et al (ref. 83) applied these concepts 
using the meteoroid mass distribution and 
the probability theory given above together 
with cratering mechanics of Gault (ref. 28). 
The number of impacts per unit area (e.g.. 
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fig. 21) were converted into “depth exca- 
vated” because each crater diameter may be 
associated with a given crater depth. Results 
of such calculations are illustrated in figure 
22. The absolute timescale is based on the 
Gault et al. micrometeoroid flux (ref. 34), 
assumed to be constant over geological times. 
Though these absolute rates of regolith turn- 
over are considered realistic for about the 
past 10* to 10® years, they are certainly not 
valid for periods > 10® years. Gault et al. 
(ref. 83) therefore also calculated the same 
data for a time variable flux; these data are 
shown in figure 23. 

The principal result of figures 21 and 23 
is of course the high turnover rate of the 
regolith surface itself, e.g., figure 22: while 
it takes approximately 10^ years to com- 
pletely turn over an 8-mm-deep zone at least 
once, the uppermost mm of the very same 
area has been turned over already 25(!) 
times; or when 99 percent of an 8-mm-deep 
layer is turned over at least once, 50 percent 
of the same surface will have already been 
turned over to a 1.4-cm depth. As a conse- 
quence, there exists a very thin surface zone, 
approximately 1 mm in thickness, in which 
extreme mixing and homogenization of com- 
ponents occurs. However, the lunar regolith 
becomes relatively quiescent rather quickly 
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Figure 22 . — The detailed turnover history of various 
regolith depths, as a function of absolute time and 
a constant flux. 


Fig:ure 23 . — Same as figure 2$, but with use of a 
meteoroid model flux that increases with geologic 
time to match the observed crater densities at the 
Apollo 12 landing site. 
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with depth; e.g., even with a meteoroid flux 
that accounts for an increase in bombard- 
ment (fig. 23) in early lunar history (=2 to 
3.8 X 10® years), a 1-m-thick layer is turned 
over only once with 99 percent confidence. 
This accounts for the observation of Russ et 
al. (ref. 78) that a major section of the 
Apollo 15 deep drill core had been residing 
completely undisturbed on the lunar surface 
for the past 500 m.y. We therefore conclude 
that due to the mass-frequency distribution 
of interplanetary matter that is vastly domi- 
nated by relative small particles in the 10-® 
to 10-^ g mass range, only an upper mm is 
thoroughly mixed before an adjacent larger 
impact event covers the area and effectively 
removes the mixing layer from the active 
reworking zone. It is thus possible to preserve 
the observed small-scale stratigraphy in the 
regolith. 

However, though each surface layer un- 
doubtedly has its peculiar surface history, it 
is not correct to conclude that each layer was 
deposited at the eventual site of recovery by 
one discrete impact event. Gault et al. (ref. 
86) and Stoffler et al. (ref. 87) demonstrated 
that the ejecta blankets of experimental im- 
pact craters in layered quartzsand targets 
have part of the original target-stratigraphy 
preserved, though in reversed sequence, i.e., 
overturned. Similar observations are also 
made around large-scale nuclear and chem- 
ical explosion craters as well as terrestrial 
impact craters (e.g., the 25-km-diameter 
Reis-structure, Germany (ref. 88). As a con- 
sequence, each regolith crater on the Moon 
will preserve — though certainly in a some- 
what degraded fashion — the original strati- 
graphic section. Therefore a variety of 
discrete layers may be excavated and rede- 
posited at the site of recovery by a large 
single impact regardless of whether they had 
drastically different exposure histories before 
this last depositional episode. 

Furthermore, processes other than direct 
deposition of impact ejecta blankets may also 
cause an apparent layering in the recovered 
regolith cores. For example: small-scale 

slumping on the walls of regolith craters may 
be a significant process. It can also be en- 


visioned that soft soil breccias ejected by a 
larger event completely desintegrate upon 
landing at significant distances from the 
primary crater. Rocks that survived such a 
landing at the end of a ballistic trajectory 
are subject to micrometeoroid erosion and 
their erosion products may be foreign to the 
new environment, thus causing a local “het- 
erogeneity” and therefore a “layer” in the 
regolith stratigraphy. Virtually nothing is 
known about the lateral dimensions of the 
regolith “layers” and it is possible that their 
areal extent is rather limited. Beyond any 
doubt, however, caution is necessary to postu- 
late that each observed layer was last de- 
posited by one discrete impact event; such 
interpretations may be grossly in error. 

LUNAR ROCK EROSION 

Studies of the grain size distribution of 
individual cratering experiments (refs. 31 
and 89) revealed that the ejecta of one given 
event are significantly more coarse-grained 
than grain sizes reported from the lunar rego- 
lith (e.g., refs. 90 and 91). Thus larger rego- 
lith components must be broken up, i.e, 
“eroded,” by small-scale cratering events. 
The visual inspection of lunar rocks both on 
lunar surface photographs and in the labora- 
tory reveals that micrometeoroid impact 
causes erosion and eventual destruction of 
rock specimens exposed to space. The micro- 
meteoroid complex operates on two different 
scales and accordingly results in two signif- 
icantly different effects, i.e., “single particle 
abrasion” and “catastrophic rupture” (refs. 
34, 59, 63, 79, and others) . 

“Single particle abrasion” is caused by 
relatively small craters, as compared with 
the overall size of a specific rock, and it re- 
sults in an effect similar to sandblasting. It 
is largely responsible for gradual mass wast- 
ing associated with a general rounding of the 
rocks (fig. 24). In contrast, “catastrophic 
rupture” is accomplished only by craters of 
relatively large size with respect to a given 
rock mass, i.e., only by impacts of sufficient 
energy to generate penetrative fracture sys- 
tems (fig. 25). 
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Horz et al. (ref. 16) simulated the “single 
particle abrasion” process via Monte Carlo- 
based computer models ; up to 10* craters 152 
to 25 000 ^ in spall diameter were produced 
on a 25-cm- surface area. Figure 26 displays 
some computer-generated profiles after a 
variety of crater numbers produced. Figure 
27 illustrates the average erosion depth as a 
function of time. Note the influence of a few, 
though large, events in particular in figure 
27, but also in figure 26. Applying a best 
estimate for the absolute micrometeoroid flux 
averaged over the past 10® years, Horz et al. 
(ref. 85) arrive at erosion rates for crystal- 
line lunar rocks of .3 to .6 mm per 10® years. 
The erosion rate for breccias may be higher, 
because of less compressive target strength 
(ref. 34). 

Figure 28 illustrates an additional result 
of the above Monte Carlo simulation relating 


to the “representative” nature of finite size 
rock chips available in the laboratory to de- 
lineate lunar surface processes. The com- 
puter iterated over the entire test surface and 
searched for the least (= shallowest) and 
most eroded (= deepest) “unit areas” that 
were defined as 5, 2, 1, .64, and .16 cm-. The 
“extremes” in erosional state are compared 
with the average of the entire area in figure 
27. The deviation from the average is a di- 
rect measure of how typical or atypical small 
lunar rock chips may be with respect to their 
parent rock. The deviations observed are 
considerable and constitute ample evidence 
that the random nature of the impact pro- 
cess has to be seriously considered in the 
analysis of discrete, finite-size rock chips. 
Unless it is demonstrated otherwise, that 
such a sample is truly “representative” of 
the parent rock, the results obtained may 



Figure 24 . — Typical lunar rock (H310) illustrating the effects of single particle abrasion. Xs indicated by the 
soil line, parts of this rock were buried in the lunar regolith. The buried portion is characterized by sharp, 
angular fracture surfaces. In contrast, the surfaces exposed to the micrometeoroid bombardment are 
abraded and significantly rounded. 
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Figure 25 . — Lunar rock 7S155 that has suffered an exceptionally large impact almost capable of catastrophi- 
cally rupturing the entire hand specimen. 


only be used with caution to delineate “aver- 
ages,” e.g., solar flare particle track densities 
to determine the absolute exposure age. 

Gault et al. (ref. 34) treated the destruc- 
tion of lunar rocks due to “catastrophic rup- 
ture.” The catastrophic breakup of rocks 
may be accomplished either by a single im- 
pact event of sufficient energy or by the cu- 
mulative effects of a number of smaller im- 
pacts; the rupture energy (B«) is cumulative 
(ref. 92). The energy required to rupture a 
rock (^ sperical body) of radius r can be de- 


scribed as 

Eg = 2.5 X 10« Sc (7) 

where Sc is the unconfined compressive 
strength of the rock in kilobars; and Eb is 
the unit energy required per gram, rather 
than total mass, of a rock of radius r. It thus 
follows that relatively less energy is re- 
quired to destroy progressively larger rock 
specimens. Figure 29 compares actual mea- 
surements of the very largest pit craters ob- 
served on lunar rocks and the relations 
expressed by equation (7). The agreement is 
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Figure 26. — Computer-generated erosion profiles of a lunar rock. According to each number of total craters 
produced, profiles taken at three different localities (Yn, Yit, Y%i) are illustrated (white — volume eroded; 
stippled — remaining rock; the vertical exaggeration is 17 times). 


good (ref. 26). Figure 29 illustrates the 
mean survival time before catastrophic rup- 
ture occurs for various hypothetical rock ma- 
terial, considering compressive strength and 
rock mass as the main variables (ref. 34). 

When the results of “single particle abra- 
sion” and “catastrophic breakup” are com- 
bined, the following conclusions emerge: 
while, for example, a 1-kg rock will sur- 
vive catastrophic desintegration for about 
3 X 10® years, it has suffered in the mean- 
time “single particle abrasion” that effec- 
tively removed a surface layer of only about 
1 to 2 mm in thickness. (See figure 30) . Thus 
“catastrophic rupture” must be considered 


the vastly superior process in obliterating 
lunar rocks. “Single particle abrasion” plays 
a minor role only, but it is still an order of 
magnitude more effective than sputtering 
processes caused by high energetic radiation 
(ref. 79) 1973. 

Conclusions 

The authors hope their efforts have demon- 
strated that the study of lunar microcraters 
has significantly contributed to our present 
understanding of the micrometeoroid com- 
plex: 
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Figure 27 . — Average erosion depth resulting from a 
Monte Carlo computer simulation. The best esti- 
mate for erosion is based on 10‘ craters, i.e., run 3. 
Notice the influence of some few — however, very 
large — craters. 

1. Contrary to popular astronomical hy- 
potheses, the micrometeoroids have 
densities from 2 to 4 g/cm®. They are 
also equant if not spherical in shape; 
forms like needles, whiskers, platelets, 
rods, etc., may safely be excluded. 

2. The mass-frequencies from 10 to 
10"^ g are in agreement with previous 
meteoroid data. However, particle 
masses as small as 10-*® g are respon- 
sible for the formation of microcraters 

< .1 pm in diameter. This result negates 
the existence of the celebrated “radia- 
tion pressure cutoff” at particle masses 

< 10 *2 g, 

3. The average micrometeoroid flux for 
the past 10® years could be established 
within a factor of 5. In agreement with 
satellite measurements, it is likely that 
the present micrometeoroid activity is 
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Figure 28 . — Extremes in deviation of erosional state 
of various, absolute surface areas (5, 2, 1., .64, and 
.16 cm') compared with the average of a 25 cm' 
surface. 

about an order of magnitude higher 
than this long-term average. 

4. Though absolute flux data do not exist 
at the moment, there is ample evidence 
that the micrometeoroid complex 
existed throughout geological time. 

5. The potential of the “lunar microme- 
teoroid detector” is not fully exhausted 
at the moment. 

The micrometeoroid complex and larger 
meteoroids are primarily responsible for the 
evolution and physical-chemical makeup of 
the lunar regolith; they effectively control 
the overall regolith growth as well as small- 
scale stratigraphy. The regolith has to be 
envisioned as a complex sequence of ejecta 
blankets that have not necessarily lost their 
integrity. The mixing, “gardening,” and 
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Figure 29. — Spall and pit diameters required for ca- 
tastrophic rupture of a given rock mass, according 
to calculations based on experimental and observa- 
tional results. Diameter, Dp (destructive) is con- 
sidered an upper limit for pit diameters observable 
on lunar rocks; Do is an experimental limit refer- 
ring to the crater diameter, i.e., spall diameter 
(D,). Ratios for D,/Dp in lunar rocks are typically 
from 3.8 to i.6. The agreement of observations on 
lunar rocks with those on experimental rupture is 
excellent. 
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Figrure 30. — Calculated mean residence time before 
destruction by catastrophic rupture for spherical 
rocks of radius r and compressive strengths (Sc), 
exposed to the micrometeoroid bombardment. 
Masses of the largest particles (Mpma) contribut- 
ing to the rupture process are indicated. 


homogenization are largely confined to the 
uppermost layer of approximately 1-mm 
thickness. Lunar rocks are effectively de- 
stroyed by micrometeoroids, with the “cat- 
astrophic rupture” process dominating the 
“single particle abrasion.” These results will 
aid not only in the interpretation of lunar 
materials but also in that of other planetary 
surfaces as well. 
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